clear;clc;
A=[9,5;2,6];
x=[1;2];
[P,D]=eig(A)
syms n
P*D^n*P^-1
ans*x